REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection 
of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

3.  DATES  COVERED  (From  -  To) 

17-10-2014 

Conference  Proceeding 

- 

4.  TITLE  AND  SUBTITLE 

How  Do  Cracks  Initiate  and  Grow  in  a  Thin  Glass  Plate?  A 
Peridynamic  Analysis 


5a.  CONTRACT  NUMBER 

W91  INF- 10- 1-0431 


5b.  GRANT  NUMBER 


6.  AUTHORS 

Florin  Bobaru,  Yenan  Wang,  Jian  Yu,  Jared  Wright,  Chian-Fong  Yen 


5c.  PROGRAM  ELEMENT  NUMBER 
622618 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 
University  of  Nebraska 

Board  of  Regents,  Univ  of  Nebraska,  Univ  of  Nebraska 
3 12  N  14th  Street 

Lincoln,  NE  68588  -0430 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS 
(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITORS  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITORS  REPORT 
NUMBER(S) 

58450-EG.27 


Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

In  this  report,  we  present  peridynamic  simulation  results  and  experimental  evidence  for  the  origin  and  evolution  of 
cracks  in  a  thin  glass  plate  that  is  impacted  near  its  center,  at  moderate  speeds  (150m/s),  by  a  small  steel  projectile. 
A  polycarbonate  thin  backing  plate  is  used  to  preserve  location  of  glass  fragments.  In  the  computational  model  the 
tape  used  on  the  sides  to  hold  the  two  plates  together  is  absent.  Upon  impact,  a  series  of  circular  cracks  form  on  the 
strike  face,  followed  by  radial  cracks  advancing  from  the  impact  site  towards  the  boundaries,  driven  by  elastic 


. ^  u a  r>  ~ ,  n . 


4-1 - u 


15.  SUBJECT  TERMS 

dynamic  fracture;  fragmentation;  brittle  fracture;  glass  plate;  multilayered  system;  transparent  protective  system;  impact 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

uu 

UU 

uu 

UU 

Florin  Bobaru 


19b.  TELEPHONE  NUMBER 
402-472-8348 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


Report  Title 

How  Do  Cracks  Initiate  and  Grow  in  a  Thin  Glass  Plate?  A  Peridynamic  Analysis 

ABSTRACT 

In  this  report,  we  present  peridynamic  simulation  results  and  experimental  evidence  for  the  origin  and  evolution  of 
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polycarbonate  thin  backing  plate  is  used  to  preserve  location  of  glass  fragments.  In  the  computational  model  the  tape 
used  on  the  sides  to  hold  the  two  plates  together  is  absent.  Upon  impact,  a  series  of  circular  cracks  form  on  the  strike 
face,  followed  by  radial  cracks  advancing  from  the  impact  site  towards  the  boundaries,  driven  by  elastic  deformations 
induced  behind  the  Rayleigh  wave.  At  the  same  time,  a  through-thickness  Hertz-crack  grows,  as  well  as  cracks 
parallel  to  the  sides  of  the  plates,  induced  by  waves  reflected  from  the  boundaries  reinforcing  surface  waves 
emanated  from  the  impact  region.  We  present  evidence  gathered  from  our  simulations  for  the  origin  and  evolution  of 
these  cracks,  and  confirm  these  results  with  fractography  experiments  of  post-mortem  samples.  The  results  provide 
evidence  of  the  predictive  capabilities  of  the  peridynamic  model  for  simulating  complex  dynamic  fracture  behavior  in 
brittle  targets  and  offer  arguments  for  what  the  sufficient  ingredients  in  a  model  for  dynamic  brittle  fracture  should 
be. 
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Plates  suspended  to  mimic 
"no  boundary  conditions". 

Impact  speeds  of  up  to 
150m/s.  No  damage 
observed  in  PC  plate. 

Tape  is  used  along 
boundaries  to  recover 
fragments  after  impact. 

Impact  location  is  off-center 

cm  closer  to  the  right). 


Setup 


Experiments  by 
J.  Yu,  C.F.  Yen  (ARL) 


0.3  cm 


•  Pictures  of  strike  face 
(taken  by  Y.  Wang  at  UNL) 

•  Tape  is  applied  on  the 
impact  face  only  post¬ 
mortem,  to  preserve 
fragments. 

•  At  higher  impact  speeds, 
fragments  are  dislodged 
and  slide  out  of  place  being 
caught  between  the  two 
plates. 


Hu,  Wang,  Yu,  Yen,  and  Bobaru, 
International  Journal  of  Impact 
Engineering  (2013) 


lental  Results 


t 

Impact  speed 
100  m/s 


Impact  speed 
150  m/s 


Strike  and  back  faces  of  glass  at  150m/s 


•  Between  33-35  major  fragments  at  150m/s  impact  speed;  many  through¬ 
thickness  cracks  are  tilted. 


•  Major  radial  cracks,  some  "branch"  before  reaching  the  boundaries. 

•  Impact  cone:  small  region  of  comminuted  material  on  strike  face,  more 
damage  on  the  back  face  of  the  plate. 

•  Major  circumferential  cracks,  and  some  surprising,  very  fine,  wispy 
lines/cracks  up  to  3.8-4  cm  diameter  around  the  impact  center. 


•  Some  through-thickness  cracks  are  parallel  to  the  sample  boundaries. 


Strike  face 


Back  face 


Questions  to  be  answered 


>  Can  we  understand  how  and  why  each  type  of  crack  system  forms? 

>  Crack  surface  fractography  can  give  indication  of  crack  speed,  crack 
growth  direction,  order  or  crack  generation 

>  Can  we  compute  such  complex  behavior? 

>  "Complexity-free  fracture/fragmentation/damage".  Can  a  single 
model,  with  minimal  inputs,  solve  the  problem  of  one  crack  and 
that  of  many  cracks  from  impact? 


The  Peridynamics  formulation 


f( 77 ,  is  a  pairwise  force  (force  density  on  particle  x  due  to  particle  x' ) 

E,  =  x'-x  (original  relative  position  ) 


A  linear  micro-elastic  material 

dw(rj  ,4) 


w(V,S) 


d  rj 

c(  Iklh* 


Ik 


+  77 


HI 


relative  elongation 


(Silling  and  Bobaru  Int  .  J .  Nonlin  .  Mech  .2005  ) 


How  cracks  form  in  peridynamics 


S.A.  Silling  (2000) 

Silling  and  Askari  (2005) 


Damage  index  number 


d  = 


n 


broken 


n 


>  d  is  between  0  and 


>  Damage  index  of  around  0.4-0. 5  localized 
along  a  line  (in  2D)  or  surface  (in  3D) 
indicates  that  a  crack  has  formed. 


Peridynamic  model 


•  No  tape  on  sides  of  the  plates  in  computations  (PC  plate  flies  away 
from  glass  after  impact). 

•  Side-taping:  rebound  of  PC  plate  causes  further  damage  in  exp. 

•  Use  15  nodes  through-thickness  of  glass;  compute  100  ps  after 
impact.  Grid  spacing  250  pm,  nonlocal  horizon  size  ~lmm. 

•  Notice  different  wave  speeds  in  glass  and  PC  (below:  strain  energy 
density  for  a  portion  cross-section  through  impact  point) 


Evolution  of  damage 


Peridynamic  results  for  impact  speed  150m/s:  damage  at 

100  ps  after  impact  _ 


•  Similar  number  of  major  fragments  as  in  experiments  (~33-35) 

•  Similar  structure  of  cone  fracture,  radial  cracks,  and 
circumferential  cracks  as  in  exp. 

•  Through-thickness  cracks  parallel  to  the  side  boundaries,  like 
observed  in  experiments. 

•  Some  cracks  "branch"  as  they  approach  the  boundaries. 

•  A  set  of  wispy  "cracks"  seen  on  the  back  face:  same  outer 
diameter  as  in  experiments  (~3.8-4  cm),  non-symmetrical. 


A  3D  view  of  Peridynamic  damage 
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strike  face,  blanking  nqt^-damaged  nodes 
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Back  face 


Tilted  view  from  the 


Strike  face 

I 


The  evolution  of  damage 

Hu,  Wang,  Yu,  Yen,  and  Bobaru,  International  Journal  of  Impact  Engineering  (2013) 


Damage  map  for  top  (top  row)  and  bottom  (bottom  row)  glass  surfaces  at 
2.95,  7.25, 15.5,  30.35,  and  99  |is  from  the  time  of  impact. 


Evolution  of  damage 


•  View  from  the  strike  face 


Two  ring-cracks,  match 
theoretical  and 


experimental  results. 
Starting  of  Hertzian  cone 


Transition  from  ring  cracks  to 
localized  radial  fractures,  that  branch, 
happens  due  to  hoop  stress  becoming 
critical  behind  the  Rayleigh  wave. 


Reflected  waves  from  the  boundaries 
generate  through-thickness  cracks 
parallel  to  the  sides 


The  evolution  of  damage  (cont'd) 


Some  cracks  start  frorntfife  boundaries 

and  join  with  radial  cracks  growing  from  the  center 


The  Hertzian  cone  transitions  to  a  crack  surface 
parallel  to  strike  face,  and  wispy  roughness 
forms  on  this  surface. 


F.P.  Bowden  and  J.E.  Field,  "The  brittle  fracture  of  solids  by  liquid  impact, 
by  solid  impact ,  and  by  shock."  Proc.  Royal  Soc  of  London,  A,  282:  331-352  (1964) 


fracture  band  radius  (in 


Peridynamic  results  for  ring  cracks 

•  For  a  3.3  mm  thick  soda-lime  glass  plate,  theory  predicts  first  radius 
of  4.2mm,  second  radius  of  6.6mm  (black  circles). 

•  PD  results:  first  radius  of  ~4mm,  second  radius  of  ~6.2mm 
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The  evolution  of  damage  (v_x  and  v_z  components) 


What  are  the  wispy  circumferential  "cracks"  ? 


Impact  speed 
150  m/s 


The  region  around  the  impact  point  (back  face) 


Damage  on 

back  side  layer  of  nodes 


See-through 
from  front 
side 


Outer  diameter  for 
wispy  "cracks":  ~3.8-4  cm 


Surprising  find  !!! 


•  Computed  wispy  "cracks"  are  not  through¬ 
thickness  cracks!  Damage  index  not  close 
to  50%  for  nodes  on  back-face. 

•  Damage  seen  as  wispy  "cracks"  is  signature 
of  roughness  on  cracks  parallel  to  strike 
face. 
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"Signature"  of  wispy  ridges  on  bottom  face 


Large  fragment  and 
glass  chip  formed 
by  the  transverse 
crack 


Wispy  roughness 
of  the  transverse  crack  that  created 


the  chip 


The  fragment  and  the  glass  chip  (view  from  bottom  side) 


The  thinner  part  of  the  chip  is  away  from  impact  center,  while  the  thicker  part  is  closer 
to  impact  center. 


Out-of-plane  nodal  velocities 


Crack  E 


Crack  G 


Boundary  cracks,  Edge-to-Center  (E2C)  cracks 


Fake  branching 
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Crease  forms  from 

Crease  initiates 

the  edge  reflection 

minor  damage 

Second  “crease”  (due  to  R-wave  reflection) 
results  in  localization  of  damage  into  cracks 


Damage  spreads  along 


edge  and  inwards 


First  “crease”  starts  a 
new  set  of  damage 
slightly  inside  the  edge 


A  3rd  set  of  E2C  cracks  is  generated  by  the  first  crease 


Symmetrical  cracks  form  on  the  lower-right  quarter  of  the  plate. 


0.3  cm 


Conclusions 


•  The  simplest  peridynamic  model  (horizon  size  is  1  mm)  with  a  relatively 
coarse  grid  (grid  spacing  is  0.25  mm)  captures  fine  roughness  of  a  crack 
surface  which  forms  when  Hertzian  fracture  cone  is  deflected  by  stress 
waves. 

•  Location  of  high  roughness  on  surface  of  crack  parallel  to  strike  face  in 
simulations  is  as  in  experiments  (shows  crack  speeds  are  computed 
correctly). 

•  Edge-to-center  cracks  produced  by  incident  waves  reinforced  by 
reflected  waves.  Fractography  confirms  one  set  of  E2C  cracks.  The  rest 
may  be  constrained  by  the  side  tape  used  in  experiments. 

•  Transition  from  ring  cracks  to  radial  cracks  (driven  by  Rayleigh  wave 
front). 

•  At  these  strain  rates/impact  speeds  and  temperatures,  nonlocal 
continuum  model  is  sufficient,  no  need  for  atomistic  or  quantum. 
Elasticity  and  damage/fracture  control  this  process. 


The  Peridynamics  formulation 


f( 77 ,  is  a  pairwise  force  (force  density  on  particle  x  due  to  particle  x' ) 

E,  =  x'-x  (original  relative  position  ) 


A  linear  micro-elastic  material 

dw(rj  ,4) 


w(V,S) 


d  rj 

c(  Iklh* 


Ik 


+  77 


HI 


relative  elongation 


(Silling  and  Bobaru  Int  .  J .  Nonlin  .  Mech  .2005  ) 


Link  to  measurable  Material  Properties 


Micromodulus  function: 

match  the  strain  energy 
for  a  homogeneous 
deformation  to  classical  elasticity 
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Bobaru  etal  (2009)  IJNME 
Convergence  to  classical ,  local  solutions 


Critical  relative  elongation:  match  the  energy/mA2  to  separate  into  two  halves 
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Silling  &  Askari,  Comp,  and  Struct.  2005  (3D) 


E.  Emmrich,  2006  (2D) 


Convergence  in  Peridynamics 
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Material  parameters 


Soda  lime  glass; 

•  Bulk  modulus:  43  GPa 

•  elastic  modulus:  72GPa 

•  density:  2440kg/mA3 

•  Poisson's  Ratio:  0.22-0.24  (0.25*) 

•  longitudinal  wave  speed:  5432m/s 

•  shear  wave  speed:  3477m/s 

•  Rayleigh  wave  speed:  3139m/s 

Bond-based  EMU  model  matches  the  bulk  modulus  43  GPa 

•  Longitudinal  wave:  5632m/s 

•  shear  wave:  3252m/s 

•  Rayleigh  wave:  2990m/s 


Dynamic  fracture/Crack  branching 


•  Experimental  observations: 

-  Branching  occurs  when  a  critical  stress  intensity  factor  is  reached  not  when  a  maximum 
velocity  is  reached. 

-  The  crack  is  at  or  near  its  maximum  velocity  for  a  considerable  path  length  before 
branching. 

-  Roughening  of  the  crack  surface  before  branching  reported  in  all  experiments 

-  Wave  propagation  can  influence  crack  branching,  curving. 


Crack  branching  in  a  single  edge  notch  specimen  of 

Bowden  et  al,  Nature  1967. 

homalite-100  (Ramulu  and  Kobayashi,  IJFM  1985). 

Possible  mechanisms  in  Crack  Branching 


•  Ravi-Chandar  and  Knauss,  IJF  (1984):  microvoids/microcracks  (glassy  polymers) 

•  Hull  (1994):  mixed-mode  loading  near  the  crack  tip  (formation  of  lances,  roughness) 


CRACK  BRANCHING  MECHANISM 


(Ramulu  and  Kobayashi,  IJFM  1985). 


/  CRACK 
$  INTERACTION 


•  Roughening  of  the  crack  surface  before  branching. 

•  Branching  occurs  when  a  critical  SIF  is  reached,  not 
when  a  maximum  velocity  is  reached. 

•  When  a  crack  reaches  a  critical  SIF  it  splits  into  two 
or  more  branches,  each  propagating  with  about  the 
same  speed  as  the  parent  crack,  but  with  a  much 
reduced  process  zone  (Ravi-Chandar,  2004). 


INSTABILITY 


Influence  of  stress  waves 


Stress  waves: 

-  influence  the  shape  and  propagation 
speed  of  dynamic  cracks. 

-  can  generate  secondary  cracks 
perpendicular  to  the  original  cracks 

-  can  induce  crack  branching  (but  are 
not  necessary  for  branching) 


Ravi-Chandar  and  Knauss,  IJF  (1984) 


Understand  crack  branching  with  peridynamics 


•  Role  of  loading  conditions  (waves). 

•  Effect  on  branching  angle,  crack  propagation 
speed,  etc. 

•  What  happens  near  the  crack  tip? 
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Sudden  loading  on  boundaries 
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9.5  jus) 


1  MPa 


2  MPa 


glass 


homalite 


Crack  propagation  speed(km/s) 


How  to  choose  the  horizon  size? 


•{>•  5=3.  Omm 

O . 5=1 ,5mm 

O  5=0. 8mm 
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Time(ps) 


Constant  speed  separation  of  boundaries 


0.2  m/s 


glass 


PMMA  (Fineberg  2006) 


Comparisons  with  experiments 


Ratio  of  speed  at  branching  to  the  Rayleigh  wave  speed. 


Loading 

Soda-lime  glass 

Homalite 

Loading  Magnitude 

V/Cr 

Loading 

Magnitude 

V/Cr 

Stresses  on  BC 

2MPa 

0.52 

IMPa 

0.46 

4MPa 

0.68 

2MPa 

0.46 

Stress  on  crack 
surface 

3MPa 

0.42 

2MPa 

0.25 
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0.44 
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0.55 

Velocity  BC 

0.06m/s 

0.58 
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0.59 
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“fish 

bones” 

Avoid  wave  interactions  with  crack  surface  loading 


0.5  MPa 


3  MPa  (branch  at  20  |lxs) 


6  MPa  (branch  at  10  jus) 


Wave  interactions  before  branching 


Wave  interactions  after  branching 


glass 


homalite 


Why/How  do  cracks  branch?  (in  this  model) 


At  higher  applied  amplitude  loading 


Migration  of  damage  away  from  the  crack  line 
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Why  do  local  models  have  trouble  with  brittle  failure? 


•  Need  to  model  the  roughness  before  branching  to  dissipate 
the  "right"  amount  of  energy  before  branching. 

•  In  dynamic  fracture,  everything  that  follows  a  "wrong" 
move  is  wrong. 

•  Models  that  use  crack  surface  tracking  (one  surface  splits  in 
two)  need  extra  conditions  to  mimic  the  actual  dissipation 
in  the  process  zone.  What  conditions? 

•  Nonlocality  helps  in  allowing  damage  to  "grow"  in  a 
preferential  direction  as  a  new  fracture  surface  with  a  much 
reduced  process  zone. 


A  "local  nonlocal"  PD  model 


Rabczuk  and  Belytschko,  IJNME  2004) 


Thickening  of  the  process  zone  before  branching 


Nonlocality  and  length-scales 


•  Peridynamic  nonlocal  damage  is  effective  in  modeling  the 
process  zone  evolution.  Once  nonlocal  scale  (horizon)  is  on 
the  scale  of  the  process  zone,  predictive  results. 

•  Migration  of  damage  away  from  the  crack  line  caused  by 
waves;  leads  to  mixed-mode  loading  conditions  and 

branching  of  the  process  zone. 

•  Dynamics  of  strain  energy  delivery  into  the  process  zone 
controls  crack  branching. 

•  Local  models  need  to  somehow  insert  this  in  the 
formulation.  Difficult  to  postulate  a  criterion,  since  it 
depends  on  local  stress/strain  conditions  (waves, 
geometry). 


